图1 黑龙江流域地理位置与海拔高度分类图
下载摘 要：温度植被干旱指数（Temperature Vegetation Dryness Index，TVDI）是指基于植被指数和地表温度相关关系估测土壤表层水分状况进而反映其所在地区的干旱程度的度量参数。本数据集覆盖的地区为黑龙江流域，时间序列为2007-2018年的4-10月（植被生长季）。该数据集是在MOD13A3植被指数（NDVI）（月度）、MOD11A2地表温度（LST）（8天合成）以及ASTER-DEM数据基础上融合温度植被干旱指数计算模型得到。数据集存储为.tif格式，空间分辨率为 1 km，由84个文件组成，数据量为362 MB（压缩为1个文件，319 MB），通过数据集与同期气象站降水和全球表层土壤水数据相关性分析表明，本产品与两者均呈负相关关系并均通过显著性检验。
关键词：黑龙江流域；生长季；温度植被干旱指数；通用NDVI-LST特征空间
DOI: https://doi.org/10.3974/geodp.2021.03.11
CSTR: https://cstr.escience.org.cn/CSTR:20146.14.2021.03.11
数据可用性声明：本文关联实体数据集已在《全球变化数据仓储电子杂志（中英文）》出版，可获取：
https://doi.org/10.3974/geodb.2021.05.02.V1 或 https://cstr.escience.org.cn/CSTR:20146.11.2021.05.02.V1.
黑龙江流域位于亚洲东北部（图1）。其南源额尔古纳河及北源石勒喀河在中国黑龙江省漠河以西的洛古河附近汇合，流域容纳俄罗斯阿穆尔州的结雅、布列亚河、阿姆贡河以及黑龙江省的松花江、乌苏里江等支流，在俄罗斯尼古拉耶夫斯克（庙街）汇入鞑靼海峡。流域面积为2,083,345.35 km2[1,2]，在世界河流流域面积排名第十[1]。作为国际上最重要的界江之一，流域跨境段长度近4,000 km，连接中蒙俄朝四国15个一级行政区，故此其内部资源环境的变化对所属的中蒙俄经济走廊沿线地区绿色发展影响重大[3–5]。近年来流域内各地区干旱事件频发，对沿线地区各领域尤其是粮食安全问题风险增大。通过多源数据及采用科学研究手段对黑龙江流域沿线旱情开展长期监测，将有助于为本流域防旱减灾对策实施及未来风险评价预测等研究工作提供科学参考及基础支持，有助于强化流域覆盖国针对该流域内各地区干旱化灾害事件的协同防治和国际合作。
图1 黑龙江流域地理位置与海拔高度分类图
下载干旱监测与预警近年来已逐渐成为学术界所关注和研究的热点议题。该领域的早期研究方法是采集站点监测的降水、气温、蒸散量、浅层土壤湿度等气象水文及土壤墒情数据并基于各类旱情评价指标进行统计分析进而量化研究区域的干旱程度[6]。由于站点受自然环境和经济成本等因素影响呈现监测数据体量有限且空间分布不均的特点，所以站点数据在监测旱情过程中所获得的结果质量无论从时间、空间、效率等层面上均受到影响。而现今遥感领域相关技术的成熟发展使得其与传统方法融合研究的设想成为可能，新方法可以高效、客观、经济地完成多尺度范围内的干旱化信息提取[7]。在基于遥感技术的旱情监测评价中，以植被指数和地表温度作为评价因子的相关指数较为常见[8]。其中温度植被干旱指数（Temperature Vegetation Dryness Index，TVDI）最初为Sandholt等[9]所提出，其指数模型原理为基于植被指数和地表温度两者的相关关系估测土壤表层水分状况进而反映其所在地区的干旱程度[10]。相关学者利用该指数进行了地理单元级、国家级、区域级、省市级等不同空间尺度的大量旱情监测应用研究。如曹晓明等应用TVDI指数方法对蒙古高原进行了近40年的长时序干旱监测，研究结果显示运用该指数可有效反映大尺度区域的旱情及其演变特征[11]。Cong等[12]基于TVDI指数计算结果与站点监测的10 cm土壤湿度数据进行组合分析验证表明利用该指标在中国东北地区进行全年度干旱动态监测是可行的。覃艺等[13]依据TVDI对2000年以来内蒙古地区生长季旱情特征变化进行监测与评估，并结合对应时期内的气温、降水等要素探讨时段内气候条件变化对其影响的趋势与幅度大小。然而当采用TVDI指数对大尺度范围地区进行干旱监测时，由于受到地表高程波动、经纬跨幅过大等影响，导致地表温度数据与实际存在较大偏差，进而影响对应特征空间的评价效果，降低TVDI指数反演表层土壤水分的数据精度[14]。针对这一问题，刘立文等[15]在构建多种常见植被指数-地表温度特征空间的过程中，采用研究区数字高程数据并对地表温度进行地形校正，减少了地形因素对地表能量平衡的影响，明显提升了基于TVDI指数反演大尺度区域土壤表层水分的效果。
《黑龙江流域生长季温度植被干旱指数1-km栅格数据集（2007‒2018）》[16]的名称、作者、地理区域、数据年代、时间分辨率、空间分辨率、数据集组成、数据出版与共享服务平台、数据共享政策等信息见表1。
| 条目 | 描述 |
|---|---|
| 数据集名称 | 黑龙江流域生长季温度植被干旱指数1-km栅格数据集（2007‒2018） |
| 数据集短名 | TVDI_AmurRiverBasin_2007-2018 |
| 作者信息 | 周业智，中国科学院地理科学与资源研究所，zhouyz@lreis.ac.cn 王卷乐 R-8881-2016，中国科学院地理科学与资源研究所，wangjl@igsnrr.ac.cn |
| 地理区域 | 地理范围为41°42′N-55°56′N，107°31′E-141°14′E，包括中国、俄罗斯、蒙古和朝鲜四国的15个一级行政区 |
| 数据年代 | 2007–2018 |
| 时间分辨率 | 月 |
| 空间分辨率 | 1 km |
| 数据格式 | .tif |
| 数据量 | 319 MB（压缩后） |
| 数据集组成 | 数据集共包括84个逐月温度植被干旱指数数据文件，文件名称为TVDI.YYYYMM.1_km_monthly.tif |
| 基金项目 | 中国科学院特别交流计划项目（Y9X90050Y2）；中国工程科技知识中心防灾减灾知识服务系统（CKCEST-2021-2-18） |
| 数据计算环境 | ENVI、ArcGIS |
| 出版与共享服务平台 | 全球变化科学研究数据出版系统 http://www.geodoi.ac.cn |
| 地址 | 北京市朝阳区大屯路甲11号100101，中国科学院地理科学与资源研究所 |
| 数据共享政策 | 全球变化科学研究数据出版系统的“数据”包括元数据（中英文）、通过《全球变化数据仓储电子杂志（中英文）》发表的实体数据集和通过《全球变化数据学报（中英文）》发表的数据论文。其共享政策如下：（1）“数据”以最便利的方式通过互联网系统免费向全社会开放，用户免费浏览、免费下载；（2）最终用户使用“数据”需要按照引用格式在参考文献或适当的位置标注数据来源；（3）增值服务用户或以任何形式散发和传播（包括通过计算机服务器）“数据”的用户需要与《全球变化数据学报（中英文）》编辑部签署书面协议，获得许可；（4）摘取“数据”中的部分记录创作新数据的作者需要遵循10%引用原则，即从本数据集中摘取的数据记录少于新数据集总记录量的10%，同时需要对摘取的数据记录标注数据来源[17] |
| 数据和论文检索系统 | DOI, CSTR, Crossref, DCI, CSCD, CNKI, SciEngine, WDS/ISC, GEOSS |
制作本数据集的遥感数据采用空间分辨率为1 km的MODIS月度植被指数产品MOD13A3[18]和8天地表温度产品MOD11A2[19]，时间范围为2007-2018年的4-10月（其中植被指数产品包括每年第91天至第274天的数据，地表温度产品包括每年第89天至第305天的数据），空间区域行列编号为h24v03、h25v03-04、h26v03-04以及h27v04。高程数据采用黑龙江流域海拔高度分类数据集[20]；气象数据采用空间分辨率为0.1°的全球逐旬表层土壤水栅格数据产品[21]以及按站点统计的逐月平均降水数据以及旬值土壤湿度数据[22,23]。
温度植被干旱指数（Temperature Vegetation Dryness Index, TVDI）是一种基于植被指数-地表温度特征空间的土壤水分监测方法[9]，其计算公式如下：
(1)
(2)
(3)
式中，T为地表温度产品的像元值，Tmin和Tmax分别为通用特征空间中NDVI所对应的最低和最高地表温度，a、c和b、d分别为NDVI-LST特征空间中干湿边的斜率与截距系数，指数取值范围为从0至1。
在NDVI-LST理论特征空间中，当研究区地表覆盖条件满足从裸地到植被全部覆盖状态、同时土壤含水量状况满足从完全干旱到全蓄水的条件时，指标散点图近似呈三角形或梯形[24]。而在实际情况下，因为受监测区域土地覆被类型及其所占面积等各方面因素的约束，常常不能完全满足理论特征空间中的所需条件，故此在基于单时期遥感影像构建特征空间的过程中，其干、湿边中的散点数据无法完全覆盖至理论特征空间边界。针对这一问题，由于监测区域的土地覆被以及土壤含水量情况在较长时间序列内会发生较大变化。故此基于长时序同期数据进行最值合成方法构建通用特征空间，该构建方案能够囊括地区内部不同时段内土地覆被和土壤湿度变化情况，能够近似得到理论特征空间地表覆盖及含水量的设定条件从而有效改善单一时段特征空间对干旱情况的反馈不足。
本数据的生产处理流程包括原始NDVI和LST数据预处理及时空序列重建，NDVI-LST通用特征空间构建及TVDI指数计算与评估等工作（图4）。数据预处理主要由分幅数据拼接与投影转换、波段提取及像元重采样等步骤组成；数据时空序列重建包括空间插补、地形校正、时间序列补齐与统一；通用特征空间构建部分主要包含单一空间干湿点提取与最值合成通用干湿边等工作；指数计算与评估的内容为将计算的研究区TVDI指数结果与降水和表层土壤水数据进行相关性分析来验证其监测适用性。
图4 黑龙江流域温度植被干旱指数数据集研发流程图
下载依据上述理论，相关学者已在不同区域的干旱监测研究中对其方案可行性及优越性予以了印证[24,25]，本研究选取2007‒2012年黑龙江流域内部国内省份（黑龙江、吉林、内蒙古）生长季代表月份（4、7、10月）的同期10 cm土壤湿度均值与基于单一和通用特征空间计算的指数结果进行相关性分析（图2），结果表明：单一时段特征空间的相关系数最小为0.352,5，最大为0.399,4。与之相比，同期通用特征空间的相关系数均有所提高，最小为0.397,1，最大为0.414,3，TVDI指数结果与实测值相关性更强。
图2 基于通用特征空间与单一特征空间的相关性结果对比图
下载本研究基于通用特征空间方法获得的LST最值变化情况以及拟合后的干湿边方程如表2和图3所示。
| 月份 | 干边方程 | 湿边方程 |
|---|---|---|
| 4月 | y = ‒20.234x + 52.094 | y = 12.216x ‒ 6.9245 |
| 5月 | y = 11.838x ‒ 6.7406 | y = 11.441x ‒ 2.8382 |
| 6月 | y = ‒6.7647x + 54.671 | y = 0.7518x + 12.232 |
| 7月 | y = ‒11.188x + 62.355 | y = ‒3.6056x +21.516 |
| 8月 | y = ‒20.381x + 68.587 | y = ‒3.8234x + 22.905 |
| 9月 | y = ‒6.9896x + 40.358 | y = 8.3332x ‒ 8.3254 |
| 10月 | y = ‒12.249x + 26.295 | y = 4.1039x ‒ 18.233 |
图3 黑龙江流域生长季各月通用特征空间变化图
下载黑龙江流域 2007–2018年生长季逐月温度植被干旱指数产品为单波段数据文件，数据名称格式为TVDI.YYYYMM.1_km_monthly.tif，其中TVDI表示数据内容为温度植被干旱指数结果；YYYYMM表示数据文件所处时间为YYYY年第MM月份；1_km表示数据集空间分辨率为1 km；monthly表示数据集时间尺度为月。数据结果在存储时放大10,000倍，即像元数值大小范围为0-10,000，使用本数据集时需乘以比例系数0.000,1。
根据TVDI指数计算的数据结果将2007‒2018年黑龙江流域逐月干旱情况划分为5级[21]：湿润（TVDI∈[0, 0.2]）、正常（TVDI∈(0.2, 0.4]）、轻旱（TVDI∈(0.4, 0.6]）、中旱（TVDI∈(0.6, 0.8]）和重旱（TVDI∈(0.8, 1.0]）。以此分级方案确定的黑龙江流域干旱等级分布结果如图5所示。
图5 基于TVDI指数的黑龙江流域2009、2013年生长季代表月干旱等级分布图
下载为反映黑龙江流域整体干旱年际变化，本文选定流域生长季的三个时期（初期4-5月，中期6-8月，末期9-10月），对各时期的TVDI指数均值进行趋势分析，分析结果如图6所示。
图6 黑龙江流域生长季各时期TVDI均值结果年际变化趋势
下载全球表层土壤水旬度数据集（Remote-sensing-based Surface Soil Moisture，RSSSM）是融合多种主、被动微波遥感数据反演出的结果，其数据精度与现今最优遥感土壤水数据产品Soil Moisture Active Passive（SMAP）为同一水平[26]。故此在本数据集的质量评估部分选用该数据集以及在研究区均匀分布的50座气象站点监测的降水量数据与本指数结果进行相关性分析（图7、图8）。在指数结果同前者逐像元级的相关性分析中，分析结果表明两者内部所有像元数据均呈现负相关关系，其中相关系数的绝对值均大于0.45，并且所有像点间的相关性全部通过了P<0.05的显著性检验。而在与降水量数据的相关性分析中，结果显示站点所在位置的TVDI指数与该数据均呈现负相关关系，相关系数绝对值的最大和最小值分别为0.558,1和0.129,1，所有数据的相关性分析结果均通过了P<0.01的显著性检验。
图7 研究区生长季代表月TVDI与RSSSM的相关系数与显著性检验结果图
下载
图8 研究区生长季各月份TVDI与降水的相关系数与显著性检验结果图
下载针对黑龙江流域这一气候变化敏感地区进行干旱监测对中蒙俄经济走廊绿色安全发展具有非常重要的现实意义，为此本文结合多源数据构建通用特征空间计算该地区2007‒2018年生长季TVDI结果，获取本时段内逐月旱情特征。由各月份通用特征空间中的干边结果可以看出，在NDVI介于0-0.2范围内多数月份的地表温度存在“爬升”现象，其原因可能是由于这些地区植被过于稀疏，近地表被风力等自然因素带走更多的热量从而导致温度较存在植被覆盖地区略低；在NDVI值大于0.8的范围内，6-10月特征空间内的干边LST值均出现明显下降的趋势，表明该时段内高植被覆盖度区域的蒸腾作用可有效削减地表附近的能量。反观湿边方程，4-5月和9-10月的地表温度的总体趋势随植被指数的增大而上升，而在生长季中期（6-8月）两者关系总体呈现负相关。两者存在差异的原因与植被的季节性蒸腾能力以及当季气温条件存在一定联系。而由黑龙江流域生长季各时期TVDI均值年际变化趋势可以看出，流域内各时期土壤含水量情况总体保持稳定，其中初期和末期的指数值均在0.6左右上下波动，而中期的指数统计值则趋近于0.4。该结果结合流域干旱等级分布图表明干旱情况相对严重区域主要集中在黑龙江流域中部及南部，其中，蒙古及中蒙交界部分在各月份内干旱等级均为最高。综合以上发现，本研究为探究黑龙江流域生长季干旱特征，为中蒙俄朝四国在对应跨境地区的干旱灾害联合防控机制提供了基础参考资料。
作者分工：王卷乐负责数据集开发总体设计工作；周业智、李凯负责数据采集处理及验证工作；王卷乐、周业智、李凯共同完成数据论文撰写。
利益冲突声明：本研究不存在研究者以及与公开研究成果有关的利益冲突。